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ABSTRACT 


A method for analytieally determining the stress 
distribution in a@ moment loaded ecireular semimonoeoque 
eylinder with eutout is presented. The teehnique employ- 
ed in the analysis consists of sombining the stringers 
and skin into a smaller number of equivalent stringers 
and shear membranes and, after determining the bending 
strain energy of the rings under various shear patterns, 
deriving the energies of the corresponding idealized 
structures. Analytieal results for the necessary patterns 
of three idealizations along with the method of applieation 
to specifie problems are presented in a form suitable 
for machine eomputation using Castigliano's Theorem. 

The comparison of the results of applieation to an exper- 
imental structure with experimental values is also set 


Perth's 
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CHAPTER I 
INTRODUCTION, STATEMENT OF PROBLEM, AND ORGANIZATION 
I. INTRODUCTION 


Because of the increasing demand for high performance 
rocket and missile systems, refined structural design and 
stress prediction methods have become very important. 

Most such techniques are, however, complicated and require 
the services of high grade engineering personnel for their 
use6é.e. On the other hand, the construction of experimental 
models for test is equally expensive and, if the results 
lead to redesign, such models are difficult to restrus tures 

Of particular interest are methods which lead to 
reasonable stress prediction for structures with cutouts 
such as are required for instrument and fuel system access 
ports. Since overall weight is a vital concern, necessary 
reinforcement of these areas must be efficiently done. If 
the stress predictions can be made with sufficient accuracy, 
consistency, and speed, the development time may be sig- 


nificantly reduced. 





II. THE PROBLEM 


It was the purpose of this study (1) to devise a 
standard method of analysis for the preliminary predietion 
of stresses in cylindrical semimonocoque shells with cutout 
when loaded in pure bending; (2) to present the results 
in such a manner that they may be used without excessive 
re-analysis to predict stresses in experimental structures; 
and (3) to compare the results of the analysis with a set 


of experimental results. 
III. ORGANIZATION 


Chapter II will contain a discussion of the previous 
related studies and the differences between them and this 
study. Also included will be a definition of the various 
special terms used. 

Chapter III is a discussion of the structural ideali- 
zation, the general method of analysis , and the steps of 
the structure analysis. 

Chapter IV presents a summary of the results of the 
analysis and the author's conclusions. | 

Appendix I illustrates the stepwise analysis of a 
representative structure. 

Appendix II presents the results of a comparison 


between analytical and experimental results. 
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CHAPTER II 


RELATED STUDIES AND DEFINITIONS OF TERMS 


I. RELATED STUDIES 


The method of analysis used in this thesis was sug- 
gested by a method for handling redundmwmt structures 
introduced by Professor W. Re Jensen in his Advaneed 
Structures course at Webb Institute of Naval Architeeture. 
Reference (1) outlines its application to a test eylinder. 
Professor Jensen made his personal notes covering addition- 
al work in the area fully available to this author. 

A closely related method of analysis is presented in 
reference (3). Basically it is a perturbation load tech- 
nique in which individual cancellation loads are placed 
at the boundaries of the desired cutout in an analytical 
structure whose external load system corresponds to the 
physical structure. The stress distribution of the struc- 
ture without cutout is presumed to be known and reference 
(4) supplies tables of coefficients to calculate the 
effect of each type of cancellation load throughout the 
structure (i.e. the resulting stress distribution in the 
unloaded structure in the vicinity of the applied unit 


perturbation loads), In application, then, perturbation 


% 





loads in magnitude and orientation necessary to cancel 
existing stresses in the loaded whole structure are ap- 
plied to the boundaries of the desired eutout and their 
associated stress distributions are superimposed on the 
existing pattern. The resultant stresses represent the 
distribution in the loaded cylinder with cutout. Refer- 
ence (5) shows an excellent agreement with the data of 
reference (2). 

In the author's opinion the two foregoing studies 
are the only ones sufficiently elose to the simplified 


analysis being sought to be pertinent. 


Il. DEFINITIONS OF SPECIAL TERMS USED 


Ssemimonocoque. As used herein, this term denotes a 
structure of uniform thickness with regularly spaced 
stiffeners in the form of ring frames and stringers. 

Orthogonal. As used herein, this term designates 
a property of loads or load systems. When more than one 
load is applied to a strueture simultaneously the effects 
of each may be considered independently if they are or- 
thogonal. The energy expression for the loaded structure 
will eontain no eross-produet terms between the several 


loads or load systems. 
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CHAPTER IIL 


IDEALIZATION AND METHOD OF ANALYSIS 


I. IDEALIZATION 


The types of structures with which this thesis is 


concerned (i.e. cylindrical semimonocoque) suggest, by 


their physical appearance, the kind of idealized structure 


which will be used. It will consist of rings, panels, 


and stringers which have the following characteristics: 


le 


The rings are loaded by shear flows only which 
act to produce ring bending. 

The panels are loaded by shear flows only. 
Their thickness is equal to the thickness of 
panels in the actual structure so that they 
will have the same effectiveness for carrying 
the shears as do panels in the actual structure. 
The stringers take longitudinal loads and 
longitudinal shear flows and have effective 
eross sectional areas which include the longi- 
tudinal load carrying capacity of the skin of 


the actual structure. 


Three idealizations, differing only in the number,of 


stringers involved, were chosen for analysis. In generals 





the number of stringers determines the complexity of the 
analysis. For this thesis, idealizations with six, eight, 


and twelve stringers were ehosen. 


II. METHOD OF ANALYSIS 


Having ehosen an idealization it was desired to set 


up the minimum number of independent shear flow and foree 
patterns which, by their subsequent combinations, would 


produce any possible loading pattern in the idealized 


structure. It is to be noted that, onee & ring shear 

flow pattern is fixed, the pattern of panel shear flows 
and stringer forees is also fixed. Figure III - 1 illus- 
trates this prineiple. The problem, then, reduces to 

that of determining the number of independent shear flow 
patterns which may be imposed upon one ring of the ideal- 
ization. Since shears must remain eonstant on the ring 
segments between stringers, must produce no net vertieal 
or horizontal forees, and must produee no net torque about 
an axis thru the center and perpendicular to the plane of 
the ring, the number of patterns is limited to the total 
number of shear segments minus the number of statie equili- 
brium conditions (i.e., in the six stringer ring, 6-3:3 
independent patterns). 


The applied load, pure bending, requires that the 


(6) 





Figure III - l 
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Resulting bay shear flow and stringer foree pattern 


Note that the top and bottom stringers have twice 
the load beeause shear flows act on both sides of 


these stringers. 
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line of intersection of its action plane with the ring 
plane be an axis of symmetry. If this line passes thru 
segments of the ring rather than thru stringer points the 
segments may have no shear. In the six stringer case, 
then, only one pattern of shear flow would be available, 
figure III-2. For this reason, as well as the fact that 
some work (Ref. (1) }) had already been done with symmetry 
thru stringers, it was decided to place the plane of action 
thru stringers. It should be noted that this makes the 
validity of the comparison with the selected experimental 
data (Ref. (2) ) somewhat doubtful. This point will be 
covered in Chapter IV and Appendix II. 

For convenience it was decided to designate individual 
ring shear flow patterns by the type of symmetry possessed 
about the vertical and horizontal axes. With S standing 
for Symmetry and A standing for Antisymmetry, the first 
letter of the hyphenated pair in a parenthesis designates 
the type of symmetry about the vertical axis while the 
second letter designates the type of symmetry about the 
horizontal axis. 

Before choosing the ring shear patterns, one more 
consideration was taken into account - orthogonality of 
load systems. As defined in Chapter II, it is a property 
of orthogonal load systems that the strain energy of a 


structure loaded by several such systems is equal to the 


(8) 





Figure III - 2 


Symmetry Axis 


This figure demonstrates that, if a symmetry axis 


is put thru the panels, the six stringer idealization 


is limited to one shear pattern, (S-S). No other 


(S-S) nor (S-A) pattern is possible which satisfies 


statie equilibrium. 
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sum of the strain energies of the strueture loaded by 

each separately. The coefficients of terms containing 

load cross products are equal to zero. Patterns with 
different symmetry are orthogonal, but, if a choice existed 
among patterns with the same symmetry, the ones chosen 

were those which were orthogonal. 

In summary, then, ring shear patterns must conform 

to the following conditions: 
1. Static equilikrium - necessary condition, 
2@e Vertical axis of symmetry - necessary condition. 
5. Orthogonality - desireable Reeder tons 

The actual ring shear patterns chosen were arrived 
at by trial and error. An interesting by-product of this 
process was the data for table III-1 which gives the 
distribution of ring shear patterns by symmetry type for 
each of the three idealizations. The proof of this table 
for the six stringer idealization is given in Appendix I, 
but it was not proven for the eight and twelve stringer 
cases. Table III-1 is presented primarily as a point of 
interest. 

Once the ring patterns were chosen, corresponding 
half structures were set up ( the structure is symmetrical 
about its centerline plane ). For a given pattern the 
fine bending strain energy was computed and expressed as 


one-half the square of the shear flow multiplied by a 


(10) 
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"flexibility factor." This factor is the same for any 
ring loaded with the given pattern. The half-structure 
was considered on a bay by bay basis as demonstrated in 
Appendix I, and the loadings on the successive bays 
starting at the half-bay (the center bay of the overall 
structure) were determined. The strain energy from stringers, 
panels, and rings was determined and the sum gave the 
strain energy of the half-structure. 

If loads are orthogonal, the sum of the strain energies 
for each pattern is equal to the strain energy of the 
idealized structure loaded by the combined patterns. Or- 


thogonality was checked for patterns of like symmetry. 


III. SUMMARY 


The steps of the analysis are: 

1. Choose an idealization. 

2. Select ring shear patterns. 

3. For each pattern set up the half-structure. 

4. For each pattern calculate the ring "Flexibility 
Factor." 

5. For eaeh pattern, on a bay by bay basis, calcu- 
late the energy matrix factors. 

Note that the combination of the patterns is not done in 


the analysis. This step is taken when the results are 


(12) 





applied to a given physical structure since only then 
will the boundary conditions be known (i.e. the required 
relationships between the pattern shears and forees and 
the forces applied to the structure). The complete analy- 
sis of a particular idealization consists of the matrices 
which correspond to the ring shear patterns for the ideal- 
ization. These are found in Chapter IV for each of the 
three idealizations. 


(13) 





CHAPTER IV 
SUMMARY OF RESULTS AND CONCLUSIONS 
Li. SUMMARY OF RESULTS 


The nature or the results indicates that the best 
method of presentation is in the form of figures and 
energy matrices. They are so presented on pages 17 to 
32 . It is intended by this form of presentation to make 
the results available in the most compact and useable 
form possible. The figures which accompany the results 
are intended to assist the reader in visualizing the 
structure to which the matrix applies. The six stringer 
idealization is presented in its entirety for both synm- 
metry patterns while, for the eight and ee stringer 
structures, only a sample bay (the half-bay, bay number 
(1) ) is illustrated for each pattern. Note that stringer 
force values are diminished as they propogate thru the 
structure and terminate in a force designated by a eapital 
letter. The end of the structure is considered as built 
into a rigid wall so that the moment of inertia of the 
end "ring" is considered to approach infinity (i.e. its 
flexibility factor is zero). The meaning of the physical 


constants is sunmarized in table IV = l. 





8 


The idealizations of this thesis are all for a nine 
bay structure because the most readily available test 
data (Ref. (2) ) were taken from such a strueture. If 
application is desired to a structure having a different 
number of bays, the matrices must be altered. The step 
by step analysis of one pattern of the six stringer idea- 
lization is presented in Appendix I so that a model is 


available for anyone wishing to make such alterations. 
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TABLE IV - 1 


TABLA OF MATRIX SYMBOL DEFINITIONS 


Shear flow. 

Radius to the neutral axis of the skin. 

Radius to the neutral axis of the skin-ring con- 
bination. 

Length of one bay., 

Thickness of the skin. 

Cross-seetional area of the idealized stringer. 
Elastic Modulus of the structure material. 

Shear Modulus of the structure material. 

Moment of Inertia of a cross-section of the skin- 


ring combination about its neutral axis, 


Underlined terms are taken from the physical struc- 


A is computed as follows: 


A=N (Ast)+1 27Rt, 
n n 


i) 


N=Number of physical stringers, 


n=Number of idealized stringers, 


Ast = Cross-sectional area of the physical stringer. 


(16) 
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In equation form, the (S-S) pattern strain energy is 
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Bach of the terms in the above energy expression 
contains a ecoeffieient times a second degree load product. 
Sueh an arrangement suggests the use of the less "eluttered" 
matrix form and, since modern digital computers are often 
adjusted to solve matrix inputs, this form was adopted. 
Conversion is aceomplished by laying out load rows and 
columns as on page (19C) and plaeing eaeh ecoeffieient in 
the bloek corresponding to the intersection of the eolumn 
of one load with the row of the other, Sinee a matrix is 
symmetric about its major diagonal (i.e. load squared 


terms), non-squared load produst positions oecur in pairs 


(19A) 








and one-half of the eoeffielent is plaeed in eaeh. For 

the same reason, only one half of the matrix need be 
filled in sinee corresponding bloeks in the other haif 

are identieal. Additional information on general matrix 
use may be found in referenee (6) while its speeifie appli- 
eation to struetural analysis is well illustrated in refer- 


enee (7). 
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IL. CONCLUSIONS 


It was originally hoped that this thesis would pro- 
vide the basis for a standardized method of preliminary 
stress analysis for the type of struseture with which it 
is coneerned. Disregarding, for the moment, the compari- 
son plots of Appendix II, the application and advantages 
of the method will first be discussed. 

Even though the simplicity of the idealization pre- 
cludes great aeeuracy, this analysis can be used in pre- 
liminary design to provide a relatively quick check on 
the adequacy of the strueture. When the design becomes 
finalized, a detailed check, using one of the more aecur- 
ate methods of stress analysis (o.g. that presented in 
reference (3)) or using an experimental model, may be made. 

In the exaet form contained herein, the analysis 
should be applied only to a nine bay strueture with load 
axis thru stringers and cutout in the center or center 
plus symmetrically placed adjacent bays. A struetural 
analyst would be required to set up the boundary equations 
defined by the physical strueture and to derive the elimi- 
nation matrix, but the calculations of matrix faetors for 
computer use and transcription of results can be aecom- 
plished by less highly trained personnel. The physical 


parameters (1.e. those items underlined on table IV - 1, 
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page (16) may be varied without rederiving the elimina- 
tion matrix so that a cheek on redesign may be rapidly 
made. 

The method of analysis may also be used to obtain 
results corresponding to struetures which vary in number 
of bays from that assumed in the thesis. Furthermore, 
by "tagging" parameters from a given bay, an analysis 
may be made of a strueture in which physical parameters 
vary from bay to bay, such as changes in stringer cross 
sections. Finally non-cutout cases may be handled by 
combining the idealizations properly. 

The major objection to the above which may be cited 
is the evidenee of poor agreement with experimental results 
seen in Appendix II. In spite of these somewhat disap- 
pointing results, it is felt that the obvious advantages 
warrant a more thorough check than time permitted. It 
would be especially interesting to have experimental stress 
data in the 180 degree area where the twelve stringer 
{idealization predicts a minimum of compression (the rever- 
sal to slight tension actually computed for one of the 
bays is considered highly unlikely). 

In summary, it is felt that, even with its relatively 
poor agreement with the only experimental data available, 
the analysis presented herein could be profitably used in 


preliminary design to quickly obtain approximate structural 
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seantlings. Further, the possible advantages seem to 


the author to warrant an experimental model set up speeli- 


fieally to verify the analytical data and an investigation 


of the effect of minor changes in analysis. 
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APPENDIX I 


DETAILED ANALYSIS OF SIX STRINGER (S-S) PATTERN 


I. AVAILABLE RING PATTERNS 


Figure AI - 1 illustrates the only three possible 


independent ring shear patterns for this idealization. 


The total number of such patterns is 


fixed by the three 


static stability conditions at six minus three equals 


three. An (A-S) pattern is not possible because a hori- 


zontal symmetry axis requires zero shear on the side 


panels while a vertical antisynmetry 
continuity of shear across the upper 
panels. With these conditions it is 
satisfy both 

ya) ieee. 
and 

2M>0. 

Sinee the load pattern requires 

axis, only the (S-S) and (S-A) shear 


able. 
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II. (S-S) STRUCTURE ANALYSIS 


Figure AI - 2 and its aecompanying calculations ill- 
ustrate the procedure for finding the bay energy as the 
sum of the energies of the elements (i.e. ring, panels, 
and stringers). 

Figures AI - 3 thru AI -7 with tables AI - 1 thru 
AI - 5 illustrate the procedure for finding the success- 
ive bay energles. Table AI = 6 represents the total 
structure energy and is the sum in each matrix block, .of 

‘ 


the corresponding blocks in the matrices of each bay. 
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Figure AI = 2 
SAMPLE BAY ENERGY CALCULATION 
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Table AI - 3 
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APPENDIX 11 
COMPARISON WITH EXPERIMENTAL RESULTS 
L. COMPARISON STRUCTURE AND PROCEDURE 


Reference (3) presents the stress data obtained by 
loading, in pure bending, an experimental cylindrical 
semimonocoque shell in which a centrally located cutout 
of varying circumferential length had been made. These 
data were chosen as the most readily available check on 
the validity of the practical application of this analy- 
sis. Aceordingly, the physical constants and combination 
equations necessary to make each analytical icealization 
conform to the experimental structure were derived. The 
eombined matrices for each idealization were solved for 
redundants on an IBM 650 computor. From the redundant 
solutions, the stringer stress values at several stations 
along the cylinder were computed. For each such station 
& plot was made of analytical and experimental stress 
values versus angular position. The zero angle position 
was vertically downward, figure ail -1, 

Because ring and skin shear values in the analytical 
structure are constant between stringer positions and 


because the angular extent of such values is so mich 
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larger than in the experimental structure, it was felt 
that plots of these quantities would not be particularly 
informative. An exception to this is the eight stringer 
case where, because the idealization provided a cutout 
size identical to that of an experimental structure, 
panel shear comparisons were made for the bay containing 
the cutout. 

Before the comparison results are presented, three 
qualifying points should be noted: 

1. The axis load symmetry in the experimental 
structure is plaeed thru the center of panels. 
All analytieal results were obtained by placing 
this axis thru stringer points - points of 
natural shear discontinuity. 

2. In both the six and twelve stringer idealizations, 
the cutout size used was 120 degrees. The 
closest experimental data are for a1350 degree 
cutout. 

3e The moment is so applied to the strueture that 
the cutout side is in tension. The same nume- 
rieal value of moment is considered to be appli- 
ed to the idealized structure as was applied 


to the physical structure. 
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Figure AIT - 1 
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II. SIX STRINGER COMPARISON 


The application of a particular idealization to a 
moment loaded physical structure is accomplished by com- 
bining the various symmetry patterns with the moment 
induced stringer loads so that the loads aétting on eut- 
away structure are iliminated. In the six stringer case, 
with a 120 degree cutout, the superposition yields the 
following equations: 

1. (2 -2g.0%) -(2#s -29s9%) — Fn =O 

E+ 850 ~ 8s@ =O 
hm is the moment induced tension load in the bottom string- 
er and equation 1. is derived from the requirement that 
this stringer have a zero load at each end of the middls 
bay. Equation 2. imposes a zero shear load on the panels 
on either side of the bottom stringer in bay l. 

The combination equations are used to eliminate 
redundants from the matrix which results from the super- 
position of the pattern matrices. These equations also 
introduce the applied moment in the form of Fm. For the 
six stringer case the final combination equations ares: 

tees = fa ~ Veale, 
Ze GsO= $20 
Figures AII - 2, 3, and 4 give the comparison plots 


with referenese data for the 150 degree cutout. 
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III, EIGHT STRINGER COMPARISON 


The three patterns must be superimposed with the 
applied moment to produce zero load in the bottom stringer 
and zero shear in the two adjacent panels. The following 


two elimination equations result. : 


i fe, = £., ~ta= ee 
Ce 4508 = 95208 - FaioX 
Figures AII - 5, 6, 7, and 8 give the sonparison 


plots with reference data for the 90 degree cutout. 


(60) 








- Aas = 4.4 “= 


Se Ean tea te ee edt ea gp lee 


« a a an 








—a 








Hy 


= 
tu 
a 
Gane 
i] 


a 
| ap 








+ 


if 
Ba 
LOSE 7 ame 
“2° tile 
TS | 
cae 


es 


+1 
ai 8 


4 . 
ict : 
oe 
ii 


pret 


4 
ond 
le 


(61 ) 





i = z 
1 Ope ae = ae ae = eee 
7 , > | 7 ] = Ae | 
: a ae 


nr 
a a ae 
| 1 ee , 





(62 


7 
‘ean 
| 2 ') 
| 


| a) | an Sc 
ae ‘eT 


= t } ; We 
i oe ee Le 
ie), -) an oa oa | 2 
° | | + 2 


: ~4orss 
| = 





a 


an 











3 


“| > wermwss DI0Mé 











| Do TESTE Se ee ee a eb 








IV. TWELVE STRINGER COMPARISON 


In order to achieve a 120 degree cutout, the three 
patterns must be superimposed with the applied moment to 
produce zero ioad in the bottom stringer and in the two . 
lowest side stringers. In addition, the four lowest 
shear panels must have zero shear. The following four 
elimination equations result. 

1 4s 1o4 = “$s2a8 - Zaiok - Ga2ok 
20 $5208 + Gs20h + Gasad 
Oe ts, = fsa + euro Ya Fm 
BPG s- 6, =\ep ta a ee 
Figures AII - 9, 10, and 1l give the scomparison 


plots with reference data for the 130 degree cutout. 
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